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ABSTRACT 


A two dimensional model of ates wines cable subjected to 
a uniform current, fixed at the sea bed and moving with a 
specified motion at the surface is developed. 

A steady state solution is developed and used as an 
initial condition for the numerical solution of the dynamic 
equations. The method of characteristics is used to carry 


out the numerical integration. 
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LI. INTRODUCTION 


The study of anchor cable dynamics in recent years has 
been of some importance to those interested in the motion 
of vessels in seaways. The mooring of buoys for navigation 
or the mooring of large tanks to hold oil or other sea- 
transported goods are further examples of the importance of 
this study. The motion of ships riding at anchor in limited 
areas would be of interest, for example, to shipping compa- 
nies that service small ports. 

In a survey [Ref. 1] done by Casarella and Parsons, it 
was mentioned that the first experimental work on mooring 
cables was probably done in 1917 by Relf and Powell. Since 
then several people have worked on the mathematical model 
of this problem. 

Through the years the model has changed to some extent. 
Some of the extensions considered have been the inclusion 
of tangential drag on the cable, the effect of surface waves, 
the effect of current, the effect of added mass, and also 
the effects of cross currents in the three dimensional 
models [Ref. 1]. Most of the studies done have been carried 
out for the static configuration; however, dynamic studies 
with arbitrary motion prescribed at the surface end of the 
cable have also been carried out. 

With the advent of fast computers, time required to get 
numerical solutions has been greatly reduced; however, 
computer time is still considerable as will be shown in this 


analysis. 
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Wilson [Ref. 2] arrived at a model very similar to the 
one used in this thesis. Wilson considered the effect of 
surface waves as well as tangential drag on the cable; 
however, the computer program used to solve the problem was 
proprietary to the company which supported his research and 
was not published. 

For the present analysis, the drag normal to the anchor 
cable only is considered and the effects of waves are 
neglected. The assumption of zero tangential drag has been 
included in the general approach to the problem and seems 
to be valid for reasonable currents and cable configurations. 
Surface waves are also neglected in this analysis as far as 
drag is concerned. The current applied is steady, uniform 
and in the horizontal plane. 

In the dynamic study, the method of characteristics is 


used to solve the equations of motion using a digital computer. 





life LHEORY 


The problem under consideration is depicted in Figure l. 
A cable of diameter d and weight W pounds per foot is anchored 
on the sea floor and is suspended from a point on the surface 
as indicated in Figure 1. A uniform current of magnitude V' 
is directed in the positive X-direction and the motion of 
the upper end of the cable may be specified in some suitable 
manner. 

Writing Newton's Law for an elemental length of cable 
as depicted in the insert of Figure 1, the following dynamic 


equation is obtained: 
T'AW = WAs" coseve= D'As' = cee (1) 
where ae Toeciec snOrmalmsacce Leration. 


Dividing equation (1) by As' and taking the limit as As'-0 


yields: 


Site ee a EC O Sia art aa ak Mae (2) 

The term M' represents the virtual mass of the cable per unit 
length and D' is defined as the cable drag per unit of length. 
T' denotes the tension in the cable and the angle Y denotes 
the angle between the tangent at a point on the cable and the 
horizontal. In summing the forces in the tangential direc- 


trom the sadaded Masis LS talsen ws zero fand as @ result the 
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actual mass is simply the cable mass. Also, the tangential 
drag is omitted as stated in the introduction. From the 
forces shown in the insert of Figure 1, the equation of 


motion for the tangential direction takes the form: 


t 
A A aR — As')cos AY - WAs'sin(¥ + =) = MA.'As* (3) 


Again, dividing through by As' and taking the limit as As'=0 


yields: 


Af ° 
ats - W sin ¥ =_ MA.' (4) 


It is necessary also to develop expressions for the 
normal and tangential accelerations in terms of the normal 
and tangential velocities. Referring to Figure 2, the normal 
and tangential acceleration components as well as the corre- 
sponding velocity components can be written in terms of the 


x and y components as 


Aes = =X cos Wee Yesin # (5) 
a = Y cos ¥ - X sin ¥ (6) 
Veg = Y sin ¥ + X cos ¥ (7) 
Vee = Y¥cos ¥ - X sin ¥ (8) 


in which x, Sey x, Y, denote the velocities and accelerations 
ie ctheaxeand ¥ directions respectively. Taking the partial 
derivatives of equations (7) and (8) with respect to time 

and using the definitions of Vis and V,' as given in equations 
(7) and (8) it is possible to show that the normal and 


tangential accelerations may be written as: 
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9Vp' a 
AW' = ee ee . 9) 
2 at’  . S§ 9! 
Hye — = - ve —_—— (10) 
Aries site. 


The primes indicate that these terms are dimensional terms. 

Dimensionless forms of these parameters will be derived later. 
Substituting equations ($) and (10) for the acceleration 

in cO anal oays (2) and (4), the following expressions for the 


two equations of motion are obtained: 








oY ao oY 
fo WwW cos Mee = aM (ee ee ey (11) 
ds! at' pen 
1 V ' 
oe--Wsiny = m (2s -v,' 2%) (12) 
ASE a Ae 





FIGURE 2 Coordinate Systems 
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PrOMemecynanlcmequattvons, Kil) and (12), 1t can bé seen 


that there are four dependent variables, V_', V.', T',¥ , all 
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of which are functions of s' and t'. In order to solve the 
equations of motion two other equations are necessary. These 
are the two kinematic relations whose derivations follow. 
With the assumption of a linear relationship between the 
stress and strain and referring to Fig. 3, it is possible to 
write the following relationship between the rate of increase 


of tension and the rate of strain, 





FIGURE 3 Tangential Forces on an Element of Cable 


As')cos A¥Y - Vg! - (Vpy't 








sin A¥] (13) 
which, upon taking the limit as As'>0, becomes: 


t a Var! 
ppot) 6 Sf See a) Ree (14) 
K ot' ao oc 
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Equation (14) states that the rate of change of tension 
equals the "Spring constant" times the rate of stretching. 
In the case of a solid cable of uniform cross section, K is 
simply AE where A is the cross sectional area and E denotes 
the modulus of elasticity. In the case of, for example, a 
chain or stranded cable the value of K would become somewhat 


more complex and would most likely be determined experimentally. 


, 8aVy aV.' 








FIGURE 4 Cable Kinematics 


A second kinematic relationship may also be introduced 
in the following manner. The angular velocity of a segment 
As' of the cable may be specified in cerns of the velocity 
of its end points. Referring to Figure 4, the angular 
velocity may be written as: 

Vy' + on as! - v,' aan 
: (15) 


q t t 1 
ot As aS As 








which becomes upon taking the limit as As'>0, 


a¥ aVy' 9 
Se nq SS (16) 
ae. as! as! 
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The four equations which have been developed, equations 
(11), (12), (14), (16), represent the complete system and 
are all that are required to determine the four dependent 


Mae Lap wese age ele, tein wuerns or the functions s? fand 


S 


t'. In summary, the basic system of equations is: 














Ti ew cosmeeep! seem (2¥n" + yon 298) (17) 
as! ane. a 
oS ew sid Yee Mies oy * 2%) (18) 
a ae ot! arte 
V | 
om oe "Ss 2 yt OMe (19) 
K ot! ds! a aa! 
qV_! 
oY = n oY 
CS = + ve) (20) 
yt! as! ys! 


The four equations which have now been specified equa- 
tions (17) through (20), may be integrated numerically to 
obtain a solution. It is possible to carry out this numerical 
integration using a finite difference formulation in the 
s'-t' plane with a rectangular grid. However, it is also 
possible, and more convenient, to use the method of charac- 
teristics as outlined by, for example, Crandall [Ref. 3]. 
Using his suggested method it can be shown that the system 
is hyperbolic and equations (17) through (20) can be trans- 
formed into the form of first order differential equations 
valid along their respective characteristics. 

According to Crandall, four other relationships in 


addition to equations (17) through (20) are needed to carry 
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out the integration by use of the method of characteristics. 


These relations are directional derivatives for the dependent 


variables Vii st V 


DV, ' 


Vee 


DY 


Dane 
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" and ¥ as defined by 





(21) 


(22) 


(23) 


(24) 


Equations (21) through (24) denote incremental changes 


associated with small displacements ds' and dt' in the 


t'-s' plane. 


In order to determine the characteristics of 


the set of equations, equations (17) through (24) are set 


in a matrix 


0 M' 0O 
0 0 0 
-1 0 0 
0 O -K 
qs rat: §6(|'0 
0 QO ds' 
0 0 0 
0 0 0 


form as follows: 


-T' M'vd 0 
0 -ViM -1 
ir 0 
ans 0 
0 0 0 
0 0 0 
ds' dt' 0 
0 QO ds' 
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© 
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iby) 


aV,'/as' 
av 1 /3at! 
Vie os. 
oWe' fat! 

o¥/os' 


owY ot" 


or'/os' 


Ql weyaore. 


-W cos ¥ - D' 


-W sin ¥ 


(25) 











Expanding the determinant of the square matrix and 


equating it to zero the following result is obtained: 


ds'4 - ( = + —)as'*at'? ao iat at'4 = 0 


Dividing by eee equation (26) becomes: 


= |X 
cs 
| 


(26) 


(27) 


Equation (27) is a quartic, the solution to which will 


produce four real characteristic directions. 


Reece 


= 
M 


2 me 
Go — 2 
i M'! 


and using these definitions equation (27) becomes: 


ase 4 = 2 Ps ase 2 2 
(GS) - (eg? +c, 1%) GED 


Solving now for ds'/dt' gives the four solutions: 


ase ' 





aes 
cca +c! 
dt! = 


18 


ZA * 
a Ce! Cy = o 


(28) 


(29) 


(30) 


(31) 


(32) 





These slopes define the four characteristics in the s‘'-t' 
plane associated with the system. The values C. and Cr 
agree with those of Wilson [Ref. 2]. The fact that there 
are four real characteristics classifies the system as 
hyperbolic. 

Having now determined the slopes of the characteristics, 
the next step is to determine the first order ordinary 
differential equations along each of these characteristics. 
To do this, the right hand side of the system matrix 
equation (Eq. 25) is substituted into any column of the 
Square matrix and this determinant is expanded and set 


equal to zero. The first determinant 1s given in eq. (33). 


Oo mt oO oO -T' M'V,' 0 -D'-W cos ¥ 
0 0 oO M 0 -M¥Vyp' -1 -W sin ¥ 


sual 0 0 0 a) rr 0 0 


= 0 (33) 


Expanding this determinant the following result is obtained: 








2 
* KC, DV," aC ae 2 DY _ eee DE" 
(K = ) ee Vers) + ( ; -Ce Neer 
C De es BR ey C 
Cc. 120 12 
= a = C.'%) W sin ¥ (34) 
C 





where Chl and C.' are defined by equations (28) and (29) 


and Cc = ds', 
ade: 


It may be noted that when C = iC,' both sides of the equation 


vanish. When C,' is substituted for C the following two 


equations, valid along the +C,' characteristic, are obtained: 


D(V,') 

(4¢,'), -——— 
Dts 
Dive) 

xe = 
Die 


Equations (35) and (36) 


Diy) 1 v(t") 
+ 
De MC,’ Dt 








D(¥) J DET) 
+ 








Dies MC, Dt? 








W sin ¥ 
(35) 
M 
-W sin ¥ 
a (36) 
M 


are the ordinary differential 


equations valid along the tC,’ characteristics; one for the 


positive slope and one for the negative slope. It is 


Piter cst in Gaeemme to heremthat the Ca slopes! are constant 


and therefore do not vary throughout the complete s'-t' 


plane. 


Equations for the C,,' slopes are determined in a similar 


manner by substituting the right hand side of the matrix 


equation system (Eq. 


(25)) into another position in the 


Square matrix; the fifth column will produce the desired 


result. After expanding, the following equation is obtained: 


D(V_,") 
Coe ea) ey e(c*-c..2) eae V.' (C,'2-C*) = 
Dt' 


Dt! 


20 


oe (@.1°-C2) (= W cos ¥ + D' 
M! 


) (37) 








Pe chtismcas> Fer e Geis er cplacedgby +C.‘|the equation 
becomes trivial. However, if +C,,' is substituted for C 
the following equations, valid along the second character- 


istics C,', are obtained: 


nN 
(+C,) Die (a=. ") dia = pe ( Wicos ? + dD ) 39 
a pt’ S nh DE Mm (38) 
Da 
gon eee 00) Sia ye A RB (SE) 
Dt' \Dhe, | M! 


Equations (35), (36), (38), and (39) represent four 
ordinary first order differential equations which are valid 


along their respective characteristics. The total deriva- 


tives RP. are defined according to eo tc ’ 
Dt! Dee mot. 2 Sk 


ds' 
where C = att denotes the slope of the characteristic in 


the s'-t' plane. Since these equations are first order 
equations they are easily integrated along the characteristic 
directions. 

With a uniform current directed as shown in Figure l, 
the drag per unit length acting normal to the cable may be 
written in terms of the fluid eget, relative EO the. capilte 


in the form 
De eesecaa(y sin ¥ + Vo“) |v" sin w+ V,"] (40) 


where op denotes the fluid density, C, denotes the drag 
coefficient, V' denotes the current magnitude and d denotes 
the cable diameter or frontal area per unit length. The 
absolute value is used in order to maintain the sign of the 


Ee latwwe sVelocity. 
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A. DIMENSIONLESS PARAMETERS 

The four prardererietic equations (35), (36), (38), and 
(39) can now be put ina dimensionless form. Various terms 
need to be defined before the actual analysis is done. 

The solution to the problem as depicted in Figure l 
depends on a number of parameters. That is, the periodic 
motion at the surface may be characterized by a period (t,) 
and a velocity amplitude (V,). Also, at steady state, the 
tension at the surface T, as well as an angle between the 
horizontal and a tangent ¥, to the cable may be specified. 
The unstretched length (L) of the cable can be used to 
normalize all length scales and the velocity amplitude may 
be used to normalize the velocities. Applying these 


definitions to the characteristic equations the following 


dimensionless parameters may be defined: 





tg (41) 
ty 
V t 

= (42) 
ve 
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(50) 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


Applying the dimensionless groupings (41) through (56) 


to equations (35), (36), (38), and (39), the characteristic 


equations can be written in the following dimensionless form: 
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(57) 


(58) 





D (V,,) DY 








(+C,) + (Vs-FCp) < = -Ccos ¥ - DIE sin ¥ + V,]|E sin ¥ + V,, | 
(39) 
D (Vy) DY . 
(-C,,) ——— + (VgtFCh) — = - C cos ¥ - DUS same, |B sive 
(60) 


Equations (57) through (60) may now be solved numerically 
on the digital computer to obtain the dynamic solution. 
However, in order to proceed with the solution it is first 
necessary to determine an initial condition as well as 
appropriate boundary Bona cions on the surface and ocean 


bottom. 


B. STEADY STATE SOLUTION 
The steady state solution may be obtained by setting Vy, 
Vs and a equal to zero in equations (57) through (60). 


Accordingly, the following steady state form of equations 


(57) and (58) is obtained: 


pol = a sin ¥ (61) 


as 


which is equivalent to equation (12) if the parameter defini- 


tions are substituted. Equations (59) and (60) become: 


FC, = C cos ¥Y + DE“sin2 y (62) 
Ss 


Equation (62) is equivalent to the equation of motion for 
the normal direction, equation (11). Since the characteris- 
tic equations agree with the equations of motion as far as 


steady state is concerned, equations (59) and (60) will be 
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the governing equations for the Steady state solution. 
Also, the fact that both sets of equations agree is an 
indication that the characteristic equations were derived 
correctly. 

Substwtuting equation (56) for C,, and dividing equation 


(61) by equation (62) provides the following expression: 





Pi 2 (63) 


This equation can be integrated analytically between the 
Pie oO le teo. and ae LOW 3 Upon integration the follow— 


ing result is obtained: 


Z 2_4 D 7 2,4 
2DE™ cos ¥-1-V1420 £ 2DE“ cos ¥ -14V1422 
(== oli) 2 sl 


_qafD?E%) * € 


| 5 
(2DE cos ¥.-1+ 144D_E re ) Ge cos ¥ -1- oe 
Cc 


C 











If there should be zero current then equation (63) would 
reduce considerably since the parameter E would drop out. 
An analysis is done in Appendix B to compare equation (64) 
without the current term to the solution obtained directly 
from equation (63) setting E equal to zero. 

Once T is known from equation (64) then the following 


expression can be used to get s: 


iC eed (65) 
AC cos ¥ + ADE’ sin? ¥ 
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Equation (65) can not be analytically evaluated, and there- 


fore it was integrated numerically on the digital computer. 


C. BOUNDARY CONDITIONS 

With the steady state solution used for the initial 
values only boundary conditions at the surface and sea bed 
need be specified. At the sea bed the normal and tangential 
velocities will be zero since the cable is anchored. With 
these two relations plus equations (58) and (60), the four 
dependent variables can be determined along the sea bed. 

At the surface, equations (57) and (59) can be used; 
however, two other equations need to be specified such that 
the four dependent variables can be calculated. These two 
relations can be specified in several ways. For example, 


the surface velocities may be written as: 


V = V 


e y Sin ¥ + V COSey (66) 


vee Vy SOK SNe sin ¥ (67) 


mie eo enc Vo may bemspecifiedc Functions Of time. 


Y 
The two relationships required at the free surface may 

also take the form of the two-dimensional equations of 

motion for a floating vessel. In this case tension would 


be related to the position, velocity and acceleration of 


the floating object. 
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III. NUMERICAL PROCEDURE 


Equations (57) through (60) represent four ordinary first 
order differential equations valid along their respective 
characteristics. In order to carry out the numerical integra- 


tion, these equations are placed in the following difference 


gor tly 
4V¥et +27. = Asin vy. Dt-v (By (68) 
Vs5 5 4G uo 1 si, mae du Gry 
Vist Veo-FCyo)¥e = (-C cos ¥,-D[E sin ¥5+ Vio I|E sin YotVi5| )Dt 
HWhot (Vo -FCyo) Yo (69) 
nst(VgstFCy3)¥, = (-C cos ¥3-DLE sin ¥,+ V3 JIE sin A Ving! Dt 
Wot (VgztFCy3) #3 (70) 
ee eee na DEV .-V + (2) TP (71) 
Pb Tw 5 Te 5 4 s4 "n4‘4 See 
ie 
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Figure 5 System of Characteristics 


Zt 





The subscripts 1,2,3,4 denote the initial point where it 
is supposed that all quantities Ver vat Te ang Ye@arne sknown. 
In order to “march out" and determine new values at the new 
location denoted by the subscript 5, the four linear equations 
(68) through (71) are solved simultaneously. Letting the 
right hand sides of equations (68) through (71) be equal to 
AA, CC, DD, and BB, respectively, the resulting expressions 


for the dependent variables at the new location 5 are: 


¥e = (CC - DD/((Vg5 - Vg3) - F(Cho + Cy3)) (72) 
2B 
Ta (A + BB +! (Vng @ Yn,))7 (=) (73) 
S 
; B 
= V  —— 
Vs5 BB + VugYs ‘on "5 (74) 
Vong = CC - (Veo - FCyo) ¥5 (75) 


The basic method of calculation indicated in equations 
(72) through (75) was used to determine values of the depen- 
dent variables Ves Vue T, and ¥ throughout the t-s plane. 

The numerical scheme utilized consisted of first dividing the 
S axis into equal subdivisions. The solution for Ves Vane Ts 
and ¥ at t=0, or initial condition, was obtained from the 
steady state solution of equations (64) and (65). The steady 
state solution provided values for Ver Moe T, and ¥ at all of 
the nodal points along the S axis. The layout in the t-s 
plane is shown in Figure 6. 


Using equations (72) through (75) a new set of values 


were then determined at the next row of points located at 
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time Dt = Ds/C. from the S axis. Once all the values of the 
four dependent variables were determined at this value of 
time, the solution was marched out another time step, Dt, so 
that values of the dependent variables were eventually 
obtained for a large area on the t-s plane. 

It may be noted in Figures 5 and 6 that the <C charac- 
teristics do not originate from a nodal point where values 
of the dependent variables are calculated. In the numerical 
procedure the method of calculation involved a linear inter- 
polation between points 1 and 4 in Figure 5 in order to 
evaluate the dependent variables at points 2 and 3. 

As indicated previously, the steady state solution was 
utilized to determine the initial values of the four dependent 
variables along the s axis. On the t axis, corresponding to 
the ocean floor, it is also necessary to apply some kind of 
boundary condition. In the example worked out, this boundary 
condition was taken as v= ‘e = 0 which corresponds to a 


faxed anchor point. 
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FIGURE 6 Layout of Characteristics 


SURFACE 
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Yee LOCUSSION OF ® RESULTS 


In order to work a numerical example, it is necessary to 
specify values for the parameters appearing in equations (41) 
through (56). For the example worked, a steel cable of uniform 
cross section with a diameter (d) of 2 inches and a modulus 
of elasticity (E) ot 30-x oe psi was selected. The density 
of the steel was 480 lbs/ft> and the cable was 1000 feet long. 
The period (t.) of the motion at the surface was 10 seconds 
and the velocity amplitude Vy was 3.0 feet/second. The 


density of seawater was taken to be 2.0 lbs-sec*/ft*. The 


Z 
Virtual mass as described by M'=M + CAP nS where Ca denotes 





the added mass coefficient [Ref. 2], and p the density of 
seawater, was.calculated to be 3.49. The added mass coeffi- 
cient used in the above calculation was taken as 1.5. The 
drag coefficient Ca was taken as 1.5 which, according to 
Wilson [Ref. 2], is valid for Reynolds numbers less than oe. 
It is now possible to calculate values for the parameters in 
equations (48) through (56). 

With an accurate steady state solution for tension having 
been calculated by equation (64), a check on the numerical 
integrations done in the characteristics method is possible. 
That is, it is possible to apply steady state boundary condi- 
tions at the free surface and recompute the steady state 
solution using the method of characteristic procedure. These 


results can then be compared with the original steady state 


solution. 


Bde 





In order to do this the steady state boundary condition at 


the surface was written in the following manner: 


i cos ween COS ve (76) 


u 
This equation states that the horizontal component of tension 
at the surface was held constant. The vertical -component was 
written as: 


S 


T sin y| 
S= a 


, Sin ie - Co, - Y_) (77) 


Equation (77) is a simplified expression representing the 
vertical component of tension as caused by a massless vessel 


floating on a still water surface in water of depth, Y. : 


O 
The dimensionless water depth Y. is defined as 1. = 
1 Oo Oo 
Lo sin¥ds|,_, and likewise the vertical span of the cable 
0 i: ] 
die time is#defined as y= f "sin ¥ ds|, . The product 
bo (0) 


Cooly, - Ys) represents a buoyant force caused by a sinking 
of the vessel a distance of (va = Ys): The constant C, is 
used to describe the interaction of the ship with the free 
surface. The "Spring rate" Cy represents the slope of the 
buoyant force-sinkage curve and is accordingly defined as 
Co = (A Yow) 2/76 where Al is the water line area, ‘ae 
is the specific weight of seawater, and L is the length of 

the cable. Using equations (76) and (77) as the surface 
boundary condition, the method of characteristics was "marched 
out" several time steps in order to compare the results 


obtained from the method of characteristics procedure with the 


initial steady state results. 





= 


The results of the above described procedure are shown 
in Figure 7. The results corresponding to t=0 were calculated 
by use of aqhasiant (64) and (65). Using these initial values, 
the method of characteristics procedure was then used to 
continue this steady state solution for some time. It may be 
seen from Figure 7 that some drift in the results took place 
due to the numerical integration by the method of characteris- 
tics when compared to the initial solution. The figure shows 
that the vessel sank slightly as the tension increased. This 
change from the initial solution is not to be construed as 
inaccuracy, however. The new steady state solution simply 
represents a second possible steady state solution. 

Since the steady state solution using the method of 
characteristics apparently provided good results, it was 
decided to add a dynamic boundary condition at the surface in 
order to induce a dynamic TESDeNaE of the cable. The boundary 
conditions at the surface, equations (76) and (77), were 
therefore modified slightly to include motion. Equation (76) 
was applied so that the horizontal component of tension was 
held constant. The vertical component of tension, however, 
was expressed in the following way: 

Sie |) ee ln Ls Col (y.-¥, Neen 2S) i) 

Figure 8 depicts the physical problem defined by the 
boundary condition equation (78). The equation describes the 
vertical component of tension as caused by a massless vessel 


Floating in a wave of elevation AL Sin 2nt measured with 
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respect to the mean water level. The term De represents the 
wave amplitude and is made dimensionless with the cable length. 
The terms within the square brackets represent the displace- 
ment of the vessel from its mean position with respect to the 
free water surface. As defined previously, Co denotes the 
slope of the buoyant force versus sinkage curve. 

Figure 9 is a graph of the dynamic response of the cable 
to the above described motion. On this graph are plotted the 
dimensionless tension and depth against a dimensionless time 
scale. Tension is made dimensionless with To which is 
defined as La WL/sin ¥) 7 te 45° being the initial steady 
state angle at the surface. The depth is made dimensionless 
with the cable length, which is 1000 feet for the example 
worked. Finally, the time is made dimensionless with the 
period of the motion at the surface which is 10 seconds. The 
To calculated from the above expression is approximately 
14,150 pounds. The results of the steady state analysis, 
Figure 7 is also shown on Figure 9. Figure 9 shows a sinu- 
soidal response which is expected from the boundary condition 
driving the motion. Although the average tension at the 
surface should be near 1.0, it can be noticed that after 
about three cycles the tension appears to be oscillating 
around a mean value which has a 6 per cent shift. This new 
mean position is equal to the one calculated from the method 
of characteristics under steady state conditions. 

Finally, it may be noted from Figure 9 that a phase shift 


exists between the position of the upper end of the cable and 


36 





((°3/4 SseTUOCTSUSUTP) AWIL 


SL°2 0S Ceeicalc 0°2 Bio Os Ge ar Ot we Gl 0S" Sou 0 
(| 26° 
Gag *e-7,  Hildad O v6° 
“Hraad SSHINOISNANIG c-7 x °x ‘Hidad SSTINOISNAWIG 
: DIWWNAG ) TLVLS AGWALS 
' + 


CORE 


yO 


; (Gene oe! 
‘NOISNaL 

4 qOWAUNS OT°T 
‘NOISNAL AOWIUNS _SSSHINOISNAWIG 
SSH INOLSNEWIG DIWYNAG 


GALVLS AGVaLS | 7 


yuzdeq pue uotsuaey otweuAq 6 FUNOLA 


ey 





the tension. Specifically, the tension tends to lag the 
displacement somewhat. This lag is caused by the dynamics 
of the phenomena; a juxtaposition of static solutions would 
necessarily show the tension and displacement to be in phase. 

The amount of computer time required to obtain numerical 
results depends on two factors, the number of subdivisions 
along the cable and the slope of the Cc. eharacteristre., “fn 
general, the number of calculations required, and consequent- 
ly, the computer time, is directly eeenont onal to the square 
of the number of subdivisions along the s- axis or the cable 
length. The size of the time step advanced in each computa- 
tion is related to the slope of the Ce characteristic 
according to Dt = Ds/C. where Ds is fixed by selection of 
the number of subdivisions of the cable and Cc. 1s defined as 
Cc. = (K/M) ? . Thus, in order to reduce the required computer 
time it is necessary to use the minimum number of subdivisions 
of the cable that will give accurate results and use small 
values of Co: Unfortunately C. 1s dependent on the cable 
physical properties and, accordingly, is fixed. 

In view of the definition of Cc. it is apparent that the 
time steps for very "stiff" cables with large elastic con- 
stants K will be small. However, for lighter cables with 
smaller elastic constants such as would correspond to nylon 
rope, the time steps will be larger and the calculation 
procedure more efficient. For the particular example worked, 
30 subdivisions on the cable were used and it took approxi- 
mately 30 minutes of computer time to compute 2.9 cycles 


corresponding to 29 seconds in real time. 
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The value of C. = (K/M) # has the physical significance 
of being the speed of propagation of a stress wave along the 
cable while the second characteristic CL = (pty)? has 
the physical significance of being the speed of propagation 
of a transverse wave along the cable. Accordingly, for all 
physically realistic cables the stress wave speed would 
normally be much larger than the transverse wave speed and 
consequently, C. would be greater than C.° In the case of 
steel cables, the stress wave speed is normally quite large 
making the time step small. Moreover, in the limiting case 
of the inextensible cable the calculation procedure breaks 
down because the time step becomes zero. However, in this 
CaceBommiucases where theycable stretching can be disregarded, 
the problem actually becomes much simpler. That is, the 
system of four equations (57-60) are replaced by the two 
equations (59) and (60) which are valid along the BC charac- 
teristics only. The *C . characteristics are lost completely 
and the system reduces to a much simpler two-characteristic 


system. 


one) : 





Vo CONGEUSIONS 


The mathematical model presented in this thesis provides 
a method for the solution of cable dynamic problems. Although 
a simplified boundary condition was employed on the surface 
for purposes of generating numerical results, the results 
obtained indicate that the method is valid and can be used in 
connection with any type of boundary condition on the surface. 

Excessive computer time is the greatest problem to be 
overcome in producing numerical results and two methods for 
reducing these were discussed. In cases where stretching can 
be disregarded the computer time can be reduced by orders of 
magnitude by utilizing the two-characteristic model since 
the time step is in this case defined by Dt = Ds/C,, and C, 
is normally small compared to Co: In cases where the cable 
is relatively stiff and yet long enough such that stretching 


cannot be disregarded, no solution to the problem of excessive 


computer time presents itself. 
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~ RECOMMENDATIONS 


The following recommendations are advanced on the basis 


of experience gained during the course of the present study: 


We 


A surface boundary condition should be applied that 
contains all the terms involved in a surface vessel's 
motion including mass, added mass, damping and buoyant 
force. 

The two-characteristic model should be developed and 
compared to the solutions using the four-characteristic 
method in order to evaluate the effect of stretching. 

A study of the effect of cable subdivision size on 
accuracy should be carried out. 

The model appearing in this thesis should be extended to 
include tangential drag, and surface wave drag on the 
cable and should be extended to the three dimensional 


Case. 
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APPENDIX B 


It is interesting to note that if the current is zero 


then equation (63) becomes: 


dr ™. ‘sin Yay 
;? = cos ¥ (78) 


This equation can be easily integrated into the following form: 
cos ¥ 
fe) 


cos’ ¥ 


By setting the current equal to zero in equation (64) it is 
not apparent that the result will be the same as equation (79). 
It 1s possible to show that equation (64) does reduce to 


equation (79) by the following method. 





Letting eS 


and using a binomial expansion for (lte)” the following 
is obtained: 


2 ol 2 


75 af oe 
(l+e)° = 1 + a 4 Te © (80) 


a 
Substituting equation (80) into the numerator of equation (64) 


produces the following numerator: 


_ 7 L “3 _ db 27m li 3 eee s a 
(- 2 5 € COS ¥. 5 € le © yt 5 e°cos ¥ + 5 
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ae ds 
The coefficient of the e«* term will produce cos ¥ since the 
other terms will be of higher order. A similar analysis can 
be done in the denominator producing a coefficient of cos ve 


Equation (64) now reduces to: 


cos ¥ 
oO 


7 ae cos !Y (82) 


which is identical to equation (79). 
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